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Maxwell  Garnett  Model  for  Dielectric  Mixtures  Containing 
Conducting  Particles  at  Optical  Frequencies 


Marina  Y.  Koledintseva,  Richard  E.  DuBroff,  and  Robert  W.  Schwartz 

University  of  Missouri-Rolla,  USA 

Abstract 

Mathematical  modeling  of  composites  made  of  a  dielectric  base  and 
randomly  oriented  metal  inclusions  is  considered.  Different  sources  of 
frequency  dependency  of  metal  conductivity  at  optical  frequencies  are  taken 
into  account.  These  include  the  skin-effect,  dimensional  (length-size) 
resonance  of  metal  particles,  and  the  Drude  model.  Also,  mean  free  path  of 
electrons  in  metals  can  be  smaller  than  the  characteristic  sizes  of 
nanoparticles,  and  this  leads  to  the  decrease  in  conductivity  of  metal 
inclusions.  These  effects  are  incorporated  in  the  Maxwell  Garnett  mixing 
formulation,  and  give  degrees  of  freedom  for  forming  desirable  optical 
frequency  characteristics  ■  of  composite  media  containing  conducting 
particles. 


I.  INTRODUCTION 

Engineering  composite  materials  with  desirable  electromagnetic  properties  for 
different  applications  from  radio  frequencies  to  the  optical  range  has  become  one  of  the 
most  important  problems  of  modern  science  and  technology. 

There  are  many  different  effective  medium  theories  used  for  modeling 
electromagnetic  properties  of  composites.  One  of  them  is  the  Maxwell  Garnett  model  [1], 
which  is  simple  and  convenient  for  modeling  due  to  its  linearity.  The  Maxwell  Garnett 
model  is  valid  for  dielectric  composites  with  dilute  conductive  phases  (below  the 
percolation  threshold).  This  is  a  model  that  implies  the  quasistatic  approximation.  Its 
main  features  are 

•  the  mixture  is  electrodynamically  isotropic; 

•  the  mixture  is  linear,  that  is,  none  of  its  constitutive  parameters  depends  on  the 
intensity  of  electromagnetic  field; 

•  the  mixture  is  non-parametric,  that  is,  its  parameters  do  not  change  in  time 
according  to  some  law  as  a  result  of  external  forces  -  electrical,  mechanical,  etc.; 

•  inclusions  are  separated  by  distances  greater  than  their  characteristic  size; 

•  the  characteristic  size  of  inclusions  is  small  compared  to  the  wavelength  in  the 
effective  medium; 

•  inclusions  are  arbitrary  randomly  oriented  ellipsoids; 

•  if  there  are  conducting  inclusions,  their  concentration  should  be  lower  than  the 
percolation  threshold. 
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Maxwell  Garnett  formalism  has  been  successfully  applied  to  engineering  microwave 
absorbing  materials  containing  carbon  particles  [2], 
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where  is  the  relative  permittivity  of  a  base  dielectric;  si  is  the  relative  permittivity  of 
the  z'-th  sort  of  inclusions;/  is  the  volume  fraction  occupied  by  the  inclusions  of  the  z'-th 
type; Njk  are  the  depolarization  factors  of  the  z'-th  type  of  inclusions,  and  the  index 

k =1,2, 3  corresponds  to  x,y,  and  z  Cartesian  coordinates.  Equation  (1)  is  generalized  for  a 
multiphase  mixture  schematically  shown  in  Figure  1.  The  constitutive  parameters  of  the 
base  material  and  inclusions  can  be  functions  of  frequency. 

Formulas  for  calculating  depolarization  factors  of  ellipsoids  and  the  table  of 
depolarization  factors  for  canonical  spheroids  (spheres,  disks,  and  cylinders)  can  be 
found  in  [3].  Equation  (1)  allows  that  within  the  same  composite  material,  particles  can 
have  different  depolarization  factors.  However,  in  reality  it  is  almost  impossible  to  have 
perfect  ellipsoidal  or  spheroidal  particles,  so,  for  any  arbitrary  shape  a  reasonable 
approximation  is  needed.  If  the  inclusions  are  thin  rods  (cylinders),  their  two 
depolarization  factors  are  close  to  Nn>2  *  1/2  ,  and  the  third  depolarization  factor  can  be 

calculated  as  in  [4],  Nn  « (l/zz)2  ln(rz).  Herein,  a  =  Hd  is  the  aspect  ratio  of  an 

inclusion,  /  is  its  length,  and  d  is  the  diameter. 

The  dielectric  properties  of  the  conducting  inclusions  are  described  by  the 
complex  relative  permittivity 


£l(jCD)  =  £’ -  j £"=£]-]- 
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the  real  part  of  which  is  much  smaller  than  the  imaginary  part  (s't  «  (J,  AT-^o))  •  In{2), 
c r,  is  the  bulk  conductivity  of  inclusions. 

It  has  been  shown  that  by  varying  conductivity,  aspect  ratio,  and  concentration  of 
inclusions  it  is  possible  to  form  the  desirable  frequency  characteristics  of  composites, 
including  shielding  effectiveness.  Application  of  the  genetic  algorithms  allows  for 
optimization  of  the  mixture  contents  to  achieve  the  highest  possible  shielding 
effectiveness  [5]. 

Mixtures  of  randomly  oriented  nanosize  conducting  particles  at  concentrations  far 
below  the  percolation  threshold  can  still  be  treated  using  Maxwell  Garnett  formalism  at 
optical  frequencies,  but  some  corrections  in  the  model  must  be  introduced.  Thus,  the 
model  should  take  into  account  peculiarities  of  the  conductive  particles  behavior  at 
optical  frequencies,  which  may  differ  considerably  from  their  behavior  at  microwave 
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frequencies.  Also,  dimensional  resonance  in  the  particles  of  inclusions  might  greatly 
influence  the  effective  permittivity  of  the  composite. 

The  structure  of  the  paper  is  the  following.  Section  II  contains  the  mathematical 
model  of  the  composite  taking  the  abovementioned  effects  into  account.  The  results  of 
calculations  based  on  this  model  are  presented  and  discussed  in  Section  III.  The 
conclusions  are  summarized  in  Section  IV. 

II.  MATHEMATICAL  MODEL 

2.1.  Renormalized  conductivity  taking  the  skin  effect  into  account 

At  the  microwave  frequencies,  the  skin  effect  in  the  conducting  inclusions  can  be 
neglected.  However,  at  the  optical  frequencies  currents  induced  by  the  electromagnetic 
waves  in  the  conducting  particles  exist  only  in  very  thin  surface  layers,  while  the  bulk  of 
the  particle  is  not  involved  in  interaction.  For  this  reason,  the  skin  effect  is  substantial  and 
influences  the  frequency  dependence  of  the  particle. 

Skin  depth  is  calculated  as 


S{co)  =  - 
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Lagarkov  and  Sarychev  [6],  according  to  [7,  Section  61],  introduced  the 
renormalized  conductivity  to  take  the  skin  effect  into  account  (see  Figure  2). 
Renormalized  conductivity  is  related  to  the  bulk  conductivity  as 


=  <?,  •  f(A) ; 


A  =  d/(2S(co))- 


m= 


1-7 

zl 
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where  J0 ,  are  the  Bessel  functions  of  the  zero  and  the  first  order,  respectively.  This 

renormalized  conductivity  appears  to  be  a  complex  value,  and  it  is  dependent  on 
frequency  through  8(a> ) .  This  is  the  basic  difference  between  the  microwave  and  optical 
behavior  of  metals:  the  conductivity  of  metals  at  optical  frequencies  is  not  constant,  it 
depends  on  frequency. 


°sk,n  (<y)  =  v'skn,  (*>)  +  jv'L  (0))  .  (5) 

The  real  part  of  conductivity  is  responsible  for  s" ,  while  the  imaginary  part  contributes 
to  s' . 
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2.2.  Conductivity  taking  into  account  small  size  of  nanorods  compared  to  mean  free 

path  of  electrons 


Since  the  inclusions  of  interest  are  nanorods,  they  are  so  thin  that  their  diameter 
could  be  smaller  than  the  electron  mean  free  path.  Then,  instead  of  the  bulk  conductivity, 
the  corrected  conductivity  should  be  taken  into  account, 

°>«  =  A(J , '  <6> 

where  the  coefficient  A  is  a  function  of  the  ratio 

A  =  f(b,ILflee),  (?) 


where  Z>,  is  the  inclusion  characteristic  size  (d  <bj  <  /  )  along  the  vector  of  the  electric 
field  acting  on  the  mixture,  and  Lfree  =  vft0  is  the  mean  free  path  for  electrons  in  the 
conductor  [8]. 

For  randomly  oriented  conducting  particles,  it  can  be  assumed  that  the  one-third 
of  all  particles  have  their  characteristic  dimension  in  the  direction  of  the  incident  electric 
field  close  to  bj  «  d ,  which  leads  to  the  decrease  of  conductivity  up  to  two  times.  One- 

third  of  all  particles  have  bt  >  Lfree  ,  so  that  their  conductivity  is  close  to  the  bulk 
conductivity.  And  the  remainder  have  an  average  characteristic  dimension  of  b,  =  ~~  > 

which  can  lead  to  some  decrease  of  conductivity.  Thus,  the  coefficient  A  may  be  even 
smaller  than  is  given  in  the  tables  in  [8],  if  surface  roughness  and  grain  size  of  metal 
inclusion  are  taken  into  account  [9-11],  For  our  computations,  we  assume  that  reasonable 
value  is  A  =  0.8 ,  that  is,  corresponding  to  a  20%  decrease  in  conductivity.  Figure  3  taken 
from  [11,  Fig.  1]  shows  that  when  the  technology  node  decreases  in  size  and  becomes 
comparable  to  the  mean  free  path  of  electrons  in  metals  (for  majority  of  metals,  such  as 
Cu,  Ag,  Pt,  Au,  Al,  etc.  it  is  in  the  range  of  10-100  nm  at  room  temperature),  the 
resistivity  of  metal  components  substantially  increases,  since  surface,  grain  boundary, 
and  barrier  layer  effects  become  substantial. 

In  such  cases,  it  is  reasonable  to  replace  the  bulk  conductivity  <x,  in  (3)  and  (4)  by 

a/ree ,  since  the  actual  conductivity  of  the  nanosize  metal  inclusion  might  be  lower  than 
the  bulk  conductivity.  Consequently, 


&M,  =crfm  • /(4);  and  5(co) 
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2.3.  Drude  model  for  metals  at  optical  frequencies 


It  is  known  that  the  frequency  dependence  of  metals  over  the  optical  frequency 
range  is  described  by  the  Drude  model  [12-15]. 


0), 


co(co-jy ) 


(9) 


where  cop  is  the  angular  plasma  frequency  for  free  electrons,  and  y  is  the  angular 
relaxation  frequency.  The  relaxation  frequency  is  y  =  l/r,  where  t  is  the  resonance 
relaxation  time.  From  (9),  real  and  imaginary  parts  of  complex  permittivity  are 


£i(jco)  =  s'i-j£"; 
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C 0 , 


2  2 
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The  imaginary  part  in  (10)  can  be  represented  as  usual 


£<  = 
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where  Drude  conductivity  is 
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(12) 


The  total  conductivity  of  a  metal  particle  is  comprised  of  “low-frequency”  conductivity 
and  Drude  conductivity  that  is  substantial  at  the  higher  optical  frequencies.  These  two 
conductivities  contribute  to  the  total  permittivity  independently,  because  their  effects  are 
separated  in  frequency. 


5 


E  (jv)  =  s',  ~ 


£i  ~  S  D  +  ^ skt 


1- 


cor 


2  2 
a)  +y 


w  skin 

COS* 


(13) 


Z,  =  Gd  +  £s*in  = 


tf)(fi>2  +  /2) 


+  - 


skin 

(OSh 


cr, 


COS, 


<J 


skin  * 


Data  for  the  plasma  frequency  and  Drude  relaxation  frequency  for  metals  can  be 
taken  from  the  graphs  and  tables  in  papers  [14,  15].  In  these  papers,  the  Drude  model 

,  .  r  .1-.  .  ri  ,  ,  ct)[rad/s]- 10”2 

parameters  co  and  y  are  given  in  [cm  J  units,  so  that  co\\l cm\  = - . 

2n-c[mls] 


2.4.  Dimensional  resonance  in  metal  particles 

Electromagnetic  resonance  in  the  particles  of  inclusions  is  another  source  of 
frequency  dependence  of  composites.  An  inclusion  particle  in  the  form  of  a  nanorod  can 
be  represented  as  an  electric  dipole  antenna.  Suppose  that  nanorods  are  non-interacting, 
corresponding  to  a  sparse  concentration  of  inclusions  in  the  base  material.  Let  us 
calculate  equivalent  conductivity  of  an  inclusion  particle.  This  conductivity  then  can  be 
used  in  Maxwell  Garnett  formulation  for  the  effective  permittivity. 

According  to  [16,17],  any  inclusion  particle  can  be  described  by  its  antenna  dipole 
input  impedance.  This  is  a  series  connection  of  radiation  resistance  Rw ,  capacitance  of 

the  wire  C,  and  its  inductance  Lj . 


Zin=Rrad  +  j®Li  +  “ 
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jcoC, 


(14) 


For  a  dipole  antenna  these  values  are  approximately  calculated  as  in  [1 8], 


Kad  =  20^  aK 


f  1  \  [Ohms];  C,  =  ,  [F];  L,=^-[ ln(2a)  - ^  | ,  [H],  (15) 

21n(2 a)  6n  (  6 
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It  is  known  that  if  the  inclusions  are  small  compared  with  wavelength,  their  input 
impedance  is  mainly  capacitive.  Radiation  resistance  determines  effectiveness  of  the 
dipole  radiation,  and  it  will  contribute  to  the  real  part  of  the  permittivity  of  the  composite. 

The  rigorous  formulation  for  radiation  resistance  and  reactance  of  a  dipole 
antenna  of  any  'length  can  be  found  in  [19,  Section  21.2],  and  these  parameters  can  be 
calculated  using  numerical  integration.  Thus,  the  input  impedance  can  be  calculated  as 
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where  the  integrand  is 
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The  distances  are  calculated  as 

R0=^dJ2f^7; 

Rt  =  V(^/2)2+(z-^)2;  (18) 

R2  =7(^/2)2+(z  +  /?)2. 


In  (1 6)-(l 8),  r/  =  l20n  j—F  (for  dielectric  composite  nef  =1),  h  =  l! 2  is  an 
'ItTc  2n  I — 

antenna  half-length,  and  k  =  —  =  —  Jse/  is  the  wave  number  in  the  effective  medium. 

Kf  K 


It  should  be  mentioned  that  sef  in  (1 5)-(  1 7)  are  taken  from  calculating  effective 

permittivity  using  the  Maxwell  Garnett  formula  (1)  before  the  dipole  antenna  effect  is 
taken  into  account. 

In  [16],  the  polarizability  of  dipole  particles  is  introduced  as 
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ee 
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(19) 


where  lef  is  the  effective  dipole  antenna  length.  For  small  dipoles,  it  can  be  assumed  that 

lef&U 2.  For  dipoles  with  /  =  / 1/2,  lef  =  0.637(2/)  =  1 .274/  [20,  Section  9.2  ]  For 

longer  dipole  antennas,  the  effective  length  is  lef  =  2.558 12 1 A .  (An  effective  length  of  a 

dipole  antenna  is  calculated  through  the  antenna  with  the  homogeneous  current 
distribution  along  the  dipole,  and  the  same  area  under  the  current  distribution  function 
along  the  dipole  length). 

The  conductivity  of  particles  due  to  the  dimensional  resonance  can  be  calculated 
by  equating  the  effective  permittivity  of  a  sparse  conglomerate  of  the  identical  conductive 
particles  (neglecting  the  base  material)  with  polarizability  ace  and  a  concentration  of  n 
particles  per  unit  volume, 
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For  cylindrical  particles, 

n  =  f,h,=2^  {Um\ 
m  l 

Then,  the  conductivity  associated  with  the  dimensional  resonance  is 

=  j(onaee  ,  [S/m].  (22) 


Substituting  (15)  and  (19)  in  (22),  one  can  get  conductivity  associated  with  the 
dimensional  resonance  as 
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(23) 


This  is  a  complex  value  ans  -  cr'res  +  jcr"res ,  and  it  can  be  added  to  the  formula 
(13),  describing  the  total  conductivity  aE.  Its  real  part  contributes  to  the  loss  in  the 
material,  and  the  imaginary  contributes  to  the  real  part  of  the  composite  permittivity, 
associated  with  phase  velocity  of  electromagnetic  waves  propagation.  This  means  that 
when  aee  is  real  (no  loss),  Z/n  is  imaginary  (purely  reactive),  then  <xrej  is  real  and 

contributes  to  loss  in  the  material.  When  aee  is  imaginary  (only  loss),  Zin  is  real  (pure 
radiation  resistance),  and  then  <jres  is  imaginary  and  contributes  to  the  real  part  of 
permittivity  of  the  composite. 


2.5.  Frequency  characteristic  of  the  base  material 

The  frequency  characteristics  of  the  base  material  also  have  an  impact  on  the 
frequency  characteristics  of  the  mixture.  Depending  on  the  base  material  used  for  making 
the  composite,  and  the  frequency  range  of  its  application,  this  base  material  can  be 
modeled  as  a  lossless  nondispersive  medium,  as  a  lossy  Debye  material,  as  a  resonance 
Lorentzian  medium,  or  as  their  superposition. 

The  Debye  frequency  characteristic  for  a  dielectric  base  material  is 
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where  esb  and  enb  are  static  and  high-frequency  (“optic”)  relative  permittivities  of  the 
material,  and  td  is  the  Debye  relaxation  time.  In  the  case  of  a  Lorentzian  material,  it  is 
described  as 
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where  co0h  is  the  angular  resonance  frequency  of  the  base  material,  and  (28 h )  is  the 
resonance  line  width  at  -3  dB  level. 

If  the  base  material  resonance  is  pronounced  and  comparatively  narrow  in  the 
frequency  range  of  interest  (narrowband  Lorentzian  material,  for  which  a>0  =  cop  >  (28) 

[22]),  it  might  have  a  great  impact  on  the  frequency-selective  properties  of  the  composite. 
If  the  dielectric  is  wideband  Lorentzian  (a>0  =  co  p  <  (28)),  then  its  behavior  is  very  close 

to  the  Debye  frequency  dependence  [22]. 

III.  RESULTS  OF  COMPUTATIONS 

A  program  in  Matlab  has  been  developed  to  calculate  reflection  and  transmission 
coefficients  from  a  slab  of  a  composite  material  containing  metal  particles  (nanorods). 
Below  there  are  the  modeled  results  that  take  into  account  all  the  abovementioned  effects. 

Figure  6  shows  the  reflection  and  transmission  coefficients  for  a  layer  of  the 
modeled  composite  material  1  pm  thick.  The  input  data  for  computations  are  the 

following:  the  bulk  conductivity  of  silver  particles  is  cr#  =  6.3  •  107  S/m;  the  coefficient  A 

taking  the  mean  free  path  into  account  is  assumed  to  be  A  =0.8;  the  mean  free  path  of 
electrons  of  silver  is  taken  as  40  nm;  the  aspect  ratio  of  the  inclusions  is  a  =  50;  the 
length  of  the  particles  /  =  1  /im;  and  the  diameter  of  the  particles  is  assumed  to  be 
d  =20  nm.  The  volume  fraction  of  the  inclusions  is/  =0.7 /a,  which  is  below  the 
percolation  threshold  estimated  as  pc=A.5/a  [6].  The  parameters  for 
polymethylmethacrylate  (PMMA)  base  material  taken  for  these  computations  are  the 
Debye  parameters:  ss  =2.2  ,  =  1.9,  and  rDehye  =  10”14  s.  It  should  be  mentioned  that 

PMMA  is  known  to  be  almost  transparent  in  the  visual  frequency  band  (98% 
transparency),  while  in  UV  range  it  almost  non-transparent  (wideband  Lorentzian 
behavior  with  coa  «  2  •  1 016  Hz).  The  dispersive  curves  for  PMMA  at  optical  and  UV 
frequencies  can  be  found  in  [21]. 

Frequency  dependence  of  reflectance,  transmittance,  and  absorbance  (magnitude, 
resonance  frequency,  and  width  of  the  resonance  line)  greatly  depend  on  the  Drude 
parameters  for  silver  -  co  and  yDrude .  Skin  effect  plays  an  important  part  for 

comparatively  low-frequency  absorption  and  frequency  selectivity.  Resonance  of 
inclusions  falls  into  the  frequency  range  of  interest  and  contributes  to  the  absorption  and 
radiation  by  inclusion  dipole  antennas.  This  contribution  is  affected  by  the  geometry  of 
inclusions. 

Figure  7  shows  the  influence  of  the  material  of  inclusions  on  the  absorbance, 
defined  as  A  =  log10(JP,  I Pmc)  in  bels.  The  Drude  parameters  for  different  metals  -  silver, 

copper,  gold,  aluminum,  and  platinum  -  are  found  in  [15]  and  presented  here  in  Table  1. 
The  base  material  was  PMMA,  and  the  inclusion  particles  had  an  aspect  ratio  of  o=50,  a 
length  of  l=\  pm,  and  a  diameter  of  <7=20  nm.  The  volumetric  fraction  of  the  inclusions 
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was/  =0.7 la.  The  composites  containing  metal  particles  with  higher  bulk  conductivity 
will  absorb  energy  more  effectively. 

Figure  8  shows  that  the  frequency  dependence  of  the  shielding  effectiveness  can 
be  also  controlled  by  varying  the  size  of  the  inclusions.  It  is  important  that  the  higher 
aspect  ratio  of  inclusions  does  not  necessarily  yield  greater  shielding  effectiveness,  as 
happens  in  the  microwave  range  [2].  At  optical  frequencies,  the  longer  particles  of 
inclusions  radiate  more  effectively  than  the  shorter  ones,  and  this  radiation  contributes  to 
the  real  part  of  the  effective  permittivity  of  the  composite,  rather  than  its  imaginary  part 
responsible  for  absorption  in  the  material. 


Table  1 

Conductivity  and  Drude  Model  Parameters  for  Some  Metals 


Metal 

Bulk  conductivity, 
10'6-  O  S/m 

10'16-  top,  1/m 

10'12-  y,  1/m 

Ag 

63.05 

1.369 

27 

Cu 

57.0 

1.122 

13 

Au 

45.55 

1.37 

40 

A1 

36.5 

2.24 

120 

Pt 

9.59 

0.78 

105 

Figure  9  demonstrates  the  effect  of  the  length  of  inclusions  on  the  absorbance. 
The  diameter  of  inclusions  for  all  the  curves  is  20  nm.  The  dimensional  resonance  in  the 
particles  corresponds  to  the  peak  of  absorption.  Thus,  for  particles  having  length  /=1000 
nm,  the  resonance  frequency  corresponds  to  a  wavelength  of  about  375  nm.  When 
particles  are  shorter,  this  peak  shifts  to  the  shorter  wavelengths.  The  absorbance 
decreases  in  the  frequency  range  of  interest  for  shorter  inclusions,  since  the  number  of  Ag 
particles  per  unit  volume  (concentration)  is  kept  constant  (4.45-10l9m'3).  Thus,  the 
volumetric  fraction  of  conducting  inclusions  decreases  when  they  are  shorter. 

Figures  10  and  11  show  the  influence  of  Ag  particles  on  the  frequency 
characteristic  of  a  mixture,  when  the  model  for  PMMA  base  material  is  nondispersive 
(Figure  10)  and  is  described  by  the  Debye  dispersion  law  (Figure  1 1).  The  base  material 
itself  is  almost  transparent  in  the  frequency  range  of  interest. 

Figure  12  contains  the  absorbance  characteristics  for  a  Lorentzian-type  base 
material,  for  a  conglomerate  comprised  solely  of  Ag  particles  of  the  same  size  and 
volumetric  fraction  as  in  Figures  10  and  11,  and  the  mixture  of  the  base  and  the  Ag 
inclusions.  The  frequency  dependence  for  a  composite  is  a  superposition  of  the  frequency 
characteristics  of  its  phases.  In  reality,  PMMA  behaves  more  like  a  Lorentzian  material  in 
the  UV  range,  where  it  becomes  nontransparent.  In  the  visual  band,  it  is  still  appropriate 
to  describe  it  as  non-dispersive  or  slightly-dispersive  with  the  Debye-like  frequency 
dependence. 
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Figure  13  contains  absorbance  characteristics  for  composites  with  different  Debye 
characteristics  of  the  base  material.  In  the  graphs  shown  in  the  figure,  the  static 
permittivity  is  a  varying  parameter.  It  is  seen  that  increase  of  pleads  not  only  to  the 
higher  absorption  in  a  wider  band,  but  also  to  the  shift  of  the  peak  of  absorption  to  the 
shorter  wavelengths.  However,  it  can  be  seen  that  substantial  increase  infiA  (from  2.2  to 

8.2)  does  not  significantly  change  the  absorption  characteristic. 

Figure  14  shows  the  frequency  dependence  of  absorbance,  when  different  models 
of  base  materials  are  taken  into  account.  There  is  almost  no  difference  in  the  curves  for 
PMMA  model  taken  as  the  Debye  (curve  1)  and  the  wideband  Lorentzian  frequency 
dependence  (curve  3).  The  size  of  the  Ag  particles  was  the  same  as  in  Figure  7  (aspect 
ratio  a= 50,  the  length  /=  1  pm,  and  the  diameter  <7=20  nm).  The  narrow  absorption  peak 
(curves  3  and  6)  can  be  obtained  using  narrowband  Lorentzian  base  material.  This  can  be 
achieved  at  co0b  =  copb  >  (2 8)  .  If  it  is  possible  to  find  a  dielectric  with  a  Lorentzian  peak 

at  visible  frequencies,  then  it  would  enhance  the  absorption  by  the  silver  particles.  This 
problem  might  be  solved,  for  example,  by  adding  Ag  particles  in  colored  optical  glass 
base  [23]. 


IV.  CONCLUSION 

Frequency-dependent  parameters  of  composites  containing  3D  randomly  oriented 
conducting  inclusions  (nanorods)  at  optical  frequencies  can  be  modeled  using  Maxwell 
Garnett  formulation,  if  the  concentration  of  inclusions  is  smaller  than  the  percolation 
threshold.  The  frequency-dependent  permittivity  of  a  mixture  mainly  depends  on  the 
conductivity  of  the  metal  inclusions.  Behavior  of  metal  inclusions  at  optical  frequencies 
greatly  differs  from  that  at  microwaves.  Nanosize  metal  inclusions  exhibit  frequency 
dependence  because  of  the  skin-effect  and  free  electron  plasma  resonance  phenomena 
described  by  the  Drude  model.  Also,  an  inclusion  particle  behaves  as  a  resonant  scatterer, 
so  its  dipole  antenna  model  should  be  taken  into  account.  In  addition,  the  mean  free  path 
of  electrons  in  metals  may  be  smaller  than  the  characteristic  size  of  metal  inclusions,  and 
this  will  substantially  decrease  the  conductivity  of  the  metals.  All  these  effects  are 
incorporated  in  the  Maxwell  Garnett  mixing  formulation,  and  give  degrees  of  freedom  for 
forming  frequency  characteristics  of  composite  media  containing  conducting  particles. 
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Figure  1.  Schematic  of  a  multiphase  mixture 
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Figure  2.  Renormalization  of  conductivity  due  to  the  skin-effect. 
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Figure  3.  Increase  of  resistivity  of  metals  (herein  -  copper)  for  elements  of  small  size 
[11]- 
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Silver;  —  ei(aj)  and  *2(60)  vs  frequency.  The  solid  line  is  the 
Drude  model.  The  data  from  Ref.  1 1  are:  Bennett  and  Bennett,  * 
for  both  —  and  *2;  Schulz,  O  for  both;  and  Hagemann  et  al X  for 

both. 


Figure  4.  Real  and  imaginary  parts  of  permittivity  for  Drude  model  of  silver  [14,  Fig  5]. 
References  are  made  on  the  paper  cited  in  [14]. 


Figure  5.  Analogy  of  an  inclusion  particle  and  a  dipole  antenna 


Reflectance  &  Transmittance 


Figure  6.  Reflectance  and  transmittance  characteristics  for  a  \  -  /u  m-thick  slab  of  a 
composite  material  PMMA-Ag. 
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Figure  7.  Absorbance  for  a  1  /u  m-  thick  slab  of  a  composite  material  -  PMMA  base  and 
different  metals  as  inclusions. 
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Figure  8.  Absorbance  control  by  varying  the  size  of  inclusions.  Volumetric  fraction  is 
/  =  0.1 1  a .  The  layer  thickness  is  1  //  m. 
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Absorbance 


Figure  9.  Absorbance  as  a  function  of  wavelength  for  a  1  /n  m-thick  composite  layer 
PMMA-Ag;  parameter  /  is  the  length  of  Ag  inclusions;  the  concentration  of  Ag  particles 

is  const  =4.45  •  1 019  m'3. 
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Absorbance 


Figure  10.  PMMA  -  nondispersive  ( ssb  =  2.2 );  Ag  particles  ( a  =  50 ;  /  =  1  //  m; 
d  =  20 nm);  concentration  nj  =  4.45  •  10“19m'3;  (/.  =  0.1  la)  for  a  \-jum  thick  slab. 
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Absorbance 


/  =  l//m;  d  =  20  nm);  concentration  nj  =  4.45-10  19 m'3  (/  =  0.7/a)  for  a  l-//m  thick 

slab. 
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Wavelength,  nm 

Figure  12.  Dielectric  -  Lorentzian  ( ssb  =  2.2  ,  =1.9,  a>0b  =  coph  =  2.05  •  1 016  rad/s 

Sb  =  1 .67 •  1 016  rad/s;)  Ag  particles  ( a  =  50 ;  l  =  ljum;  d  =  20 nm);  concentration 
ni  =  4.45  ■  1 0“19  m'3  ( /  =  0.7  /  a )  for  a  1  -  //  m  thick  slab. 


Absorbance 


Figure  13.  Ag  particles  (a  =  50;  /  =  1  ^m;  d  -  20 nm);  concentration 
ni  =  4.45  'lO-19 1/m3  (/,  =  0.7/a)  for  a  1-^m  thick  slab  for  five  different  base 
materials,  denoted  as  Debyel  through  Debye  5  corresponding  to  ssb  =  2.2 , 4.2,  8.2,  20.2, 
and  50.2.  The  parameters  sxh  =  1 .9  and  r  =  1  •  1 0"14 1/s  are  the  same  for  all  the  five  cases. 
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Absorbance 


Figure  14.  Ag  particles  (a  =  50;  /  =  1  /im;  d  =  20 nm);  concentration  nf  =  4.45  - 1 0_19m' 
3  (/  =  0.7 / o )  for  a  \- ju m  thick  slab  with  five  different  base  materials:  1-  PMMA  - 
Debye;  2-  Lorentzian  1 :  ssb  =2.2,  exh  =1.9,  coQb  =  copb  =  2.05  ■  1016  rad/s;  <5  =  1 .67  •  1016 
rad/s;  3 -Lorentzian  2:  ssb  =2.2,  smh  =1.9,  co0b  -copb  =2.05-1015  rad/s;  8b  =  1 .67  •  1 016 
rad/s;  4  -  Lorentzian  3 :  ssh  =2.2,  sxh  =1.9,  co0b  =  copb  =2.05-1 015  rad/s;  8b  =  1 .67  •  1 015 
rad/s;  5 -Lorentzian  4:  ssb  =2.2,  sv:h  =1.9,  co0b  =a>pb  =5.0-1 015  rad/s;  8b  =1.67-1015 
rad/s;  6  -  Lorentzian  5:  ssb  =2.2,  sxb  =1.9,  co0b  =a>pb  =10.0-1015  rad/s;  8b  =  1 .67  •  1 015 
rad/s. 
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